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Abstract : Our research focused on eutrophication control and species 
screening for riparian zone vegetation restoration in the upstream reach 
of the Hun River. We studied 16 hardwood plant species to investigate 
nutrient concentrations and nitrogen and phosphorus accumulations. 
After about 120 days of growth in pots, these 16 species varied in dry 
matter biomass, ranging from 15.13 to 637.16 g. Total nitrogen (TN) and 
total phosphorus (TP) concentrations and distribution in roots, stems and 
foliage differed both within and between tested species. Mean TN and TP 
accumulation ranged from 0.167 to 14.730 g per plant and from 0.016 to 
1.20 g, respectively. All 16 species, but especially Lespedeza bicolor , 
Robinia pseudoacacia and Sorbaria sorbifolia had strong potential to 
remove TN and TP from soil and could be widely utilized for the restora¬ 
tion of destroyed riparian zones in northeast China. 
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Introduction 

Eutrophication is a widespread and increasing problem in water 
resources in many countries (Bennett et al. 2001). Agricultural 
non-point source pollution remains the greatest global contribu¬ 
tor to eutrophication (Corwin et al. 1997; Dabrowski et al. 2002). 


Project funding: This work was funded by the major National Science 
and Technology project "Water Pollution Control and Management" 
(2012ZX07202008) of China and the National Science and Technology 

Pillar Pro gram (2012BAC05B05). _ 

The online version is available at http:// www.springerlink.com 

Shuai Yu *' 2 , Wei Chen Xingyuan He(^2)', Zhouli Liu 1 
Yanqing Huang 1 

1 State Key Laboratory of Forest and Soil Ecology, Institute of Applied 
Ecology, Chinese Academy of Sciences, Shenyang 110016, People’s 
Republic of China; 2 University of Chinese Academy of Sciences, Bei¬ 
jing 100049, People’s Republic of China. 

E-mail: hexv@,iae.ac.cn; oncehere88@gmail.com 

Corresponding editor: Zhu Hong 


The riparian zone receives and retains large amounts of nutrient 
inputs from farmland. These inputs can indirectly lead to eutro¬ 
phication of rivers (Giese et al. 2003; McClain et al. 2003). Ri¬ 
parian zone vegetation, through interception and storage, plays 
an important role in protecting the river water from eutrophica¬ 
tion (Lowrance et al. 1997; Hazlett et al. 2008; McBroom et al. 
2008). Nitrogen concentration in runoff water can be reduced by 
65%-100% after passing through a riparian zone forest 
(Spoelstra et al. 2010). Depending on the plant species, plant 
nutrient uptake has been shown to account for 3%-47% of ni¬ 
trogen removal and 3%-60% of phosphorus removal from runoff 
water in the riparian zone (Cooke 1992; Tanner 1996; 
Kuusemets et al. 2001). The riparian zone fonns an appropriate 
environment for nutrient removal (Lowrance et al. 1984; Lee et 
al. 2009). Previous studies have examined more than 50 aquatic 
plants (Borin and Salvato 2012) and herbs (Mcjannet et al. 1995; 
Yu et al. 2014). On an annual basis, forested riparian buffer 
strips have proven more effective at reducing nitrogen concen¬ 
trations in streams than herbaceous buffers (Lowrance et al. 1984; 
Osborne and Kovacic 1993; Hefting et al. 2005). Poplar forest in 
the riparian zone proved more effective (99% retention of N0 3 ) 
than grass (84% retention of N0 3 ~) during winter months 
(Haycock and Pinay 1993). 

Selecting suitable woody plants for restoration of polluted, 
damaged riparian zones is an effective and efficient measure for 
controlling non-point-source agricultural inputs of nutrients 
(Bedford et al. 1999). However, few studies directly compared 
the effect of various plant species under similar conditions for 
extended periods of time to identify those which most efficiently 
accumulate nitrogen and phosphorus. This is especially true in 
cold temperate regions. 

To fill the above information gap, we studied 16 native woody 
plants of 9 families to assess their performance in removing ni¬ 
trogen and phosphorus from runoff water. We quantified Total 
nitrogen and total phosphorus concentrations, and biomass ac¬ 
cumulations of the plants in different tissues. This research was 
undertaken to aid in predicting the response of species to eutro¬ 
phication and providing statistical support for plant species 
screening for the restoration of eutrophicated riparian zones. 
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Materials and methods 

Plant culture 

Based on species lists for riparian zone forests along the Hun 
River, we selected sixteen woody plant species of northeastern 
China for study: Syringci reticulate, Prunus padus, Robinia 
pseudoacacia, Pterocarya stenoptera, Juglans mandshuriea, 
Berberis dielsiana, Sambucus williamsii, Salix matsudana, 
Quercus mongolica, Rosa davurica, Euonymus alatus, Acer 
truncation, Populus alba, Lespedeza bicolor, Ulmus pumila, and 
Sorbaria sorbifolia. The uniform plants (three-year-old seedlings) 
used in this study were bought from a nursery and cultivated in 
pots (five pots per species) with 8-kg un-contaminated soil 
(Meadow burozem soil, 6.82 pH, 2.53% organic carbon, 4.35% 
organic matter, 1.88 mg-g' 1 total nitrogen, 0.24 mg-g" 1 total 
phosphorus. After one month of growth in pots, a screening ex¬ 
periment was conducted from May to October 2012 at Shenyang 
Arboretums of Chinese Academy of Sciences (41°54' N, 
123°35’ E). Plants were grown under nutrient controlled condi¬ 
tions with nutrient solution added. 

Nutrient solution was used to simulate the high nutrient input 
levels which typically occurred in agriculture runoff. The stan¬ 
dard solution contained 56 mg-L" 1 of nitrogen (as NH 4 N0 3 ) and 
62 mg-L" 1 of phosphorus (as NaH 2 P0 4 -2H 2 0), a control group 
was set up (CK: 0 mg-L" 1 of nitrogen and 0 mg-L" 1 of phospho¬ 
rus). The nutrient addition was applied once a week. The ex¬ 
periment was replicated three times. 

Nutrient analysis 

At the end of the growing season (October), the plant samples 
were harvested, washed and separated into roots, stem and foli¬ 
age, heated in an oven at 90°C for 30 min and dried at 65°C to 
constant weight. Oven-dried materials were milled and passed 
through a 100-mesh (0.149 mm) nylon sieve (Huafeng, Zhejiang, 
China), and then stored in jars prior to laboratory analyses (Lu 
1999). 

Subsamples were digested for total nitrogen (TN) and total 
phosphorus (TP) measurement according to the sulfuric ac¬ 
id-hydrogen peroxide (H 2 S0 4 -H 2 0 2 ) method (Son and Gower 
1992). TN was quantified using the semi-micro Macro Kjeldahl 
method (Ruizheng Kjeldahl nitrogen analyzer KDY-600D, 
Shanghai, China). TP was quantified using the molybdenum 
antimony-ascorbic acid colorimetric method (MADAC) (SHI- 
MADZUUV-1800 spectrophotometer, Japan). 

Statistical analyses 

Average values and standard errors (S.E.) were calculated by 
Microsoft Office Excel 2007 for all data. Statistical procedures 
used in this study were perfonned using SPSS (Version 16.0, 
SPSS Inc. 2007). Standard one-way analyses of variance 
(ANOVA) were used to test significance of differences among 
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and between the sixteen tested species and between roots, stems 
and foliage, with respect to TN, TP and biomass. Duncan multi¬ 
ple range test was employed to show the variation in TN and TP 
between species. Spearman’s correlation analysis was used to 
detennine the correlationships between tissues and nutrition. 
Significant and extremely significant differences were set as p < 
0.05 and p < 0.01, respectively. Hierarchical cluster analysis was 
used to classify plants into different groups based on the nutrient 
distribution in tissues. Origin 8.0 was used to draw figures. 

Results 

Total biomass 

At the end of the research, the total biomass per plant was com¬ 
pared between the 16 tested species (Fig. 1). After about 120 
days of growth, total biomass/plant differed significantly by 
species and treatment (p <0.05, n=3, Fig. 1). The biomass of CK 
ranged from 4.20 to 237.89 g in the order: Quercus mongolica < 
Ulmus pumila <Pterocarya stenoptera <Sambucus williamsii 
<Euonymus alatus < Berberis dielsiana <Prunus padus <Rosa 
davurica <Acer truncatum <Juglans mandshuriea < Syringa 
reticulate <Salix matsudana <Populus alba <Sorbaria sorbifo¬ 
lia < Robinia pseudoacacia <Lespedeza bicolor. The T1 (treat¬ 
ment) biomass per plant value ranged from 15.13 to 637.16 g in 
the following order: Pterocarya stenoptera < Quercus mongo¬ 
lica <Ulmus pumila <Berberis dielsiana <Sambucus williamsii 
< Prunus padus <Rosa davurica < Euonymus alatus < Acer 
truncatum<Juglans <mandshuriea < Syringa reticulate <Salix 
matsudana < Populus alba <Sorbaria sorbifolia <Lespedeza 
bicolor < Robinia pseudoacacia. T1 treatment resulted in greater 
biomass than CK. Robinia pseudoacacia and Lespedeza bicolor 
yielded more biomass than other species. 



Plant species 

Fig. 1: The CK and T1 biomass of 16 plant species. A, B, C indicate T1 
the significance different at p <0. 05; a, b, c indicate CK the significance 
different at p <0.05. Error bars denote 1 SE (n=3). Syringa reticulate (Sr), 
Prunus padus (Pp), Robinia pseudoacacia (Rp), Pterocarya stenoptera (Ps), 
Juglans mandshuriea (Jm), Berberis dielsiana (Bd), Sambucus williamsii, 
Salix matsudana (Sm), Quercus mongolica (Qm), Rosa davurica (Rd), Euo¬ 
nymus alatus (Ea), Acer truncatum (At), Populus alba (Pa), Lespedeza bicolor 
(Lb), Ulmus pumila (Up) and Sorbaria sorbifolia (Ss). 
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Variation in TN and TP concentration among tissues 

For all tested species, TN and TP concentrations varied between 
roots, stem and foliage (Fig. 2). TN (p <0.05) and TP (p <0.05) 
concentrations in foliage were significantly greater than in roots 
and stems. There was no consistent pattern in the distribution of 
N or P concentrations in stems and roots. The plant tissues in T1 
treatment had higher N and P concentrations than that did CK. 
The foliage of Sambucus williamsii and Berberis dielsiana had 
the highest TN (35.56 mg-g' 1 ) and TP (5.49 mg-g' 1 ) concentra¬ 
tions. The average TN concentration in roots, stems and foliage 
was 12.11 mg-g' 1 , 12.96 mg-g' 1 and 22.46 mg-g' 1 ; The TP con¬ 
centration in roots, stems and foliage was 1.58 mg-g' 1 , 1.51 
mg-g' 1 and 2.51 mg-g' 1 , respectively. 

Accumulation of TN and TP in plants 

According to the analysis of TN and TP accumulations in the 
plants, there were significant differences in the roots, stem, foli¬ 
age and among species. As shown in Tables 1 and 2, mean TN 
accumulated in the whole plant ranged from 166.65 to 14729.73 
mg ( Quercus mongolica <Pterocarya stenoptera <Ulmus pu- 
mila< Rosa davurica <Prunus padus < Sambucus williamsii < 
Berberis dielsiana <Acer truncatum <Euonymus alatus <Jug- 
lans mandshurieci <Syringa reticulate <Salix matsudana 
<Populus alba <Sorbaria sorbifolia < Lespedeza bicolor < 
Robinia pseudoacacia) and TP ranged from 16.06 to 1203.54 mg 
(;Quercus mongolica <Pterocatya stenoptera <Ulmus 
pumila<Euonymus alatus <Rosa davurica < Sambucus william¬ 
sii < Prunus padus <Berberis dielsiana <Acer truncatum <Salix 
matsudana < Populus alba <Juglans mandshuriea <Syringa 
reticulate <Sorbaria sorbifolia <Lespedeza bicolor <Robinia 
pseudoacacia). Species Rp showed by a wide margin the highest 
accumulated TN and TP. 


Spearman correlation analysis indicated that TN and TP ac¬ 
cumulations in different tissues were positively correlated (Table 
3). The Correlation coefficient ranged from 0.81 to 0.96. 
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Fig. 2: Total Nitrogen and Total Phosphorus concentration (dry, wt) of 
plants. Syringa reticulate (Sr), Prunus padus (Pp), Robinia pseudoacacia 
(Rp), Pterocarya stenoptera (Ps), Juglans mandshuriea (Jm), Berberis dielsi¬ 
ana (Bd), Sambucus williamsii , Salix matsudana (Sm), Quercus mongolica 
(Qm), Rosa davurica (Rd), Euonymus alatus (Ea), Acer truncatum (At), 
Populus alba (Pa), Lespedeza bicolor (Lb), Ulmus pumila (Up) and Sorbaria 
sorbifolia (Ss). 


Table 1: Comparison of the TN accumulations (mg) in different tissues (roots, stem and foliage) of the sixteen tested tree species (means ± sd) 


Total Nitrogen 

Roots (mg) 

Stem (mg) 

Foliage (mg) 

Total (mg) 

Syringa reticulata 

607.67446.65c 

438.0846.93b 

230.5443.47a 

1276.29450.12A 

Prunus padus 

227.91420.69b 

205.4140.17ab 

149.19416.92a 

582.52437.43A 

Robinia pseudoacacia 

3898.984546.31a 

4861.1941151a 

5969.5741169.87a 

14729.7341926.56C 

Pterocarya stenoptera 

141.7746.18a 

65.22416.24a 

74.24428.32a 

281.22438.39A 

Juglans mandshuriea 

508.55448.41b 

245.39420.33a 

461.4445.14b 

1215.34423.59A 

Berberis dielsiana 

256.8411.13b 

179.4418.89a 

209.61418.77ab 

645.81448.79A 

Sambucus williamsii 

293.5844.6c 

112.7540.38a 

228.16417.8b 

634.49413.58A 

Salix matsudana 

590.114209.36a 

337.42454.43a 

418.654182.68a 

1346.19481.11A 

Quercus mongolica 

68.13429.46a 

40.3143.71a 

58.2243.35a 

166.65429.82A 

Rosa davurica 

238.99414.11a 

181.34416.28a 

77.1449.21b 

497.47439.59A 

Euonymus alatus 

417.21410.66b 

202.3245.72a 

193.78414.47a 

813.3141.91A 

Acer truncatum 

255.09413.34a 

167.0449.83a 

379.43442.44b 

801.56465.6A 

Populus alba 

517.97432.11b 

814.64120.35b 

120.1242.51a 

1452.684149.95A 

Lespedeza bicolor 

998.744119.38a 

1788.724197.67b 

1466.894190.63ab 

3956.234358.43B 

Ulmus pumila 

78.6844.14a 

87.64416.89a 

176.23425.32b 

342.55446.35A 

Sorbaria sorbifolia 

643.41425.35b 

704.77438.06b 

354.05432.57a 

1702.24495.99A 


Notes: a, b, c in the same row indicate the significance different at/? < 0. 05. A, B, C in the same line indicates the significance different at/? < 0. 05. 
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Table 2: Comparison of the TP accumulations (mg) in different tissues (roots, stem and foliage) of the sixteen tested tree species (means ± sd) 


Total Phosphorus 

Roots (mg) 

Stem (mg) 

Foliage (mg) 

Total (mg) 

Syringa reticulata 

120.2849.1b 

59.2440.35a 

40.7743.14a 

220.1645.61AB 

Prunus padus 

38.3444.44b 

30.6940.72ab 

21.3242.61a 

90.246.33AB 

Robinia pseudoacacia 

256.99416.08a 

502.794123.27a 

441.384102.13a 

1203.544191.09C 

Pterocarya stenoptera 

9.76±0.56a 

641.57a 

8.2142.44a 

28.2643.45A 

Juglans mandshuriea 

60.0744.33a 

43.7342.82a 

59.147.92a 

162.9440.77AB 

Berberis dielsiana 

24.2641.71a 

30.0444.32a 

36.7344.75a 

91.32410.77AB 

Sambucus williamsii 

44.3440.56b 

21.1741.09a 

23.3941.58a 

88.9140.07AB 

Salix matsudana 

41.5340.83a 

80.7248.82b 

21.5541.33a 

143.6149.32AB 

Quercus mongolica 

7.2742.87a 

4.840.52a 

3.7340.2a 

16.0643.59A 

Rosa davurica 

39.6941.54c 

32.3241.54b 

14.640.75a 

86.5142.33AB 

Euonymus alatus 

37.441.06b 

14.7940.16a 

17.0940.61a 

69.2741.84A 

Acer truncatum 

48.842.83b 

26.6440.24a 

56.5644.31b 

132.0747.38AB 

Populus alba 

65.0942.44a 

65.7647.9a 

16.4141.13b 

147.1746.59AB 

Lespedeza bicolor 

94.41416.4a 

119.29421.75a 

115.48421.53a 

328.39421.13B 

Ulmus pumila 

15.3640.04a 

13.0841.87a 

14.7242.48a 

43.1444.39A 

Sorbaria sorbifolia 

114.2242.52c 

91.5942.21b 

43.5944.14a 

249.3543.84AB 


Notes: a, b, c in the same row indicate the significance different at/> < 0. 05. A, B, C in the same line indicates the significance different at p < 0. 05. 


Distribution of TN and TP in tissues 

The accumulated quantities of TN in roots, stems and foliage 
accounted for 23-48%, 18-56% and 8-52%, respectively, of the 
total accumulations (Fig. 3). The proportions of TP in roots, 
stems, and foliage were 21-55%, 20-56% and 11-43%, respec¬ 
tively, of the total accumulations (Fig. 3). 


Table 3: Spearman correlation coefficient of TN and TP accumulations 
in different tissues. 


Tissues 

Nutrient 

Roots 

Stem 

Foliage 

TN TP 

TN 

TP 

TN 

TP 




TN 

1 0.92** 

0.91** 

0.89** 

0.89** 

0.88** 

Roots 

TP 

1 

0.87** 

0.90** 

0.84** 

0.84** 


TN 


1 

0.95** 

0.85** 

0.81** 

Stem 

TP 



1 

0.85** 

0.84** 

Foliage 

TN 

TP 




1 

0.96** 

1 


Notes: ** Correlation is significant at the 0.01 level (2-tailed). TN is Total 
Nitrogen; TP is Total Phosphorus. 

Based on the distribution of nutrient accumulations, the tested 
species were clustered into three distinct groups by Hierarchical 
cluster analysis (Fig. 4). The first group, including only two spe¬ 
cies, had a relatively higher proportion of nutrients (more than 
50%) in stems; the second group included 9 species in which 
roots had the highest percentages of nutrients (about 45%). Foli¬ 
age and stems shared almost the equal proportions. The third 
group, including 5 species, shared nearly equal percentages of 
nutrients in roots, stems and foliage. Species of the same families 
were clustered into the same groups (except Quercus mongo- 
lica). 
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Fig. 3 : Total nitrogen and total phosphorus distribution in roots, stem and 
foliage. Syringa reticulate (Sr), Prunus padus (Pp), Robinia pseudoacacia 
(Rp), Pterocarya stenoptera (Ps), Juglans mandshuriea (Jm), Berberis dielsi- 
ana (Bd), Sambucus williamsii , Salix matsudana (Sm), Quercus mongolica 
(Qm), Rosa davurica (Rd), Euonymus alatus (Ea), Acer truncatum (At), 
Populus alba (Pa), Lespedeza bicolor (Lb), Ulmus pumila (Up) and Sorbaria 
sorbifolia (Ss). 

Discussion 

Interspecific variation 

Biomass accumulations of nutrients varied by species (Fig. 1). 
These differences were largely detennined by their physiology 
(i.e., net assimilation) and morphology (Ma et al. 2010), as many 
species have various mechanisms for adaptation, including ad- 
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justments of growth rate, modifications of plant structure (Li et 
al. 2007). Plants with high biomass accumulate more nutrients in 
their tissues (Jiang et al. 2011). Robinici pseudoacacia and Les- 
pedeza bicolor• yielded the highest biomass and accumulated the 
highest quantities of nutritients. Kyambadde et al. (2004) and 
Iamchaturapatr et al. (2007) report that species like Robinia 
pseudoacacia and Lespedeza bicolor are more suitable for ripar¬ 
ian zone restoration owing to higher N and P removal from wa¬ 
ter. 

Dendrogram using Average Linkage (Between Groups) 

Rescaled Distance cluster combine 


0 5 10 15 20 25 



Fig. 4: Euclidean distance clustering tree. Syringa reticulate (Sr), Primus 
padus (Pp), Robinia pseudoacacia (Rp), Pterocatya stenoptera (Ps), Juglans 
mandshuriea (Jm), Berberis dielsiana (Bd), Sambucus williamsii, Salix ma- 
tsudana (Sm), Quercus mongolica (Qm), Rosa davurica (Rd), Euonymus 
alatus (Ea), Acer truncation (At), Populus alba (Pa), Lespedeza bicolor (Lb), 
Ulmuspumila (Up), and Sorbaria sorbifolia (Ss). 

Boyd (1970, 1978) observed large interspecific variation in 
nutrient concentrations in aquatic plants. In his study, it is un¬ 
clear whether such differences were related to environmental 
nutrient levels or to the different absorption rates of the various 
species because the plants were collected from the field. How¬ 
ever, in our study, plants were of the same age and were culti¬ 
vated in the same environment conditions. In the absence of 
other factors such as disturbance, interspecific variation in tissue 
TN and TP was recorded in our study. This result was in accor¬ 
dance with research on 41 wetland plants reported by McJannet 
(1995). 

In our study, each species showed dramatically different TN 
and TP concentrations and accumulations between tissues (roots, 
stems and foliage). Nutritients were recorded at higher concen¬ 
trations in foliage than in roots and stems. There were, however, 
no consistent differences in nutrient concentrations between 
roots and stems. This result was consistent with the results of 
Zhu et al. (2011). Li et al. (2013) reported nutrition distribution 
as foliage> stems> roots in 30 common plant species grown in 
the hydro-fluctuation belt of Baihua Reservoir in Guizhou prov¬ 


ince, China. This was because both the structure and function 
differ by tissue type. Leaves contain photosynthetic tissues 
whose metabolism is active, while roots and stems are storage 
tissues that transport water and nutrients. Stems and roots, which 
were primarily composed of cellulose, have lower nutrition de¬ 
mand (Shan et al. 2011). 

We recorded positive correlation between TN and TP. Ni- 
inemets and Kull (2003), however, found no correlation between 
TN and TP in plant species in a wooded meadow and a bog, 
probably because concentrations were similar in all species. In 
fact, both N and P can stimulate growth or other processes be¬ 
cause TN supply affects how efficiently TP is acquired and used, 
and vice versa (Treseder and Vitousek 2001; Gusewell et al. 
2003). 

The function of riparian zone plants in nutrition removal 

TN and TP concentration of the 16 woody plants ranged from 
8.38 to 35.56 mg-g' 1 and from 0.87 to 5.49 mg-g' 1 , respectively. 
TN and TP were accumulated to quantities ranging from 111.6 to 
14729.73 mg and 16.06 to 1203.54 mg, respectively. Of the 16 
woody plants, Lespedeza bicolor, Robinia pseudoacacia, and 
Sorbaria sorbifolia had absorbed most TN and TP from soil and 
stored most in tissues. These three species are recommended as 
preferred restoration plants for the main purpose of TN and TP 
removal in the riparian zone. Of course, the above results showed 
only the absorption and storage ability of the tested plant species. 
In other words, the data in Table 2 and Table 3 were just a part of 
the total removal effect by the whole riparian ecosystem. The 
riparian plant community provides a suitable environment for TN 
and TP removal (Wu et al. 2011). Other mechanisms, such as 
rhizosphere microbial activity and physical processes, could also 
contribute to the removal of most pollutants (Brix 1987; 
Gottschall et al. 2007). Consequently, the nutrient treatment and 
removal capacities of woody plants examined this study would, 
if grown in the wild, undoubtedly far exceed the results pre¬ 
sented here. Due to the complexity of riparian ecosystems, more 
research is needed to learn more about the processes and mecha¬ 
nisms in natural situations. 
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